function ret = SIQR (vec)
    ret=(prctile(vec,75)+prctile(vec, 25))/2;
end